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1. Introductory remarks 

One-dimensional spin-i XY model in a transverse field defined by the Hamilto- 



nian 



is known as the simplest quantum many-body system for which many statistical me- 
chanical calculations can be performed exactly, i.e., without making any simplifying 
approximations. For more than forty years this model has been a standard testing- 
ground for checking various conjectures or new calculation schemes and approaches 
in statistical mechanics and condensed matter physics. The aim of the present paper 
is to elucidate some recent results derived for spin-| XY chains and to foresee some 
further problems which are attractive to study. The interest in spin-| XY chains 
may be enforced nowadays because of the progress in material sciences (see, for ex- 
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ample, the recent report on CS2C0CI4, the compound which is a good reahzation of 
the famous spin-| isotropic XY chain [E]). 



2. Dynamic properties in fermionic picture 

Spin- 1 XY chains contain a hidden symmetry which was discovered by applying 
the Jordan- Wigner transformation: the system of interacting spins can be de- 
scribed in terms of noninteracting spinless fermions (E.Lieb, T.Schultz, D.Mattis). 
As a result many (although by no means all) statistical mechanical calculations 
can be performed rigorously. As an example of notorious problems in the statis- 
tical mechanics of spin-i XY chains we may mention the analysis of the time- 
dependent correlation functions of x {y) spin components (s^(t)s^_,_^), s^(t) = 
exp(ii7t)s^exp(-ii/t), ((. . .)) = Tr(exp(-/3if)(. . .))/Tr(exp(-/3/7)). Since the rela- 
tion between the x spin component attached to a certain site and the on-site creation 
and annihilation operators of fermions is nonlocal and involves the occupation- 
number operators of fermions at all previous sites, the problem of applying the 
Wick-Bloch-de Dominicis theorem to a product of a huge number of multipliers 
arises. The result can be written compactly as the Pfaffian of an antisymmetric ma- 
trix (generally speaking of huge sizes) constructed from the elementary contractions 
and hence a further analytical analysis becomes not simple. The problem has been 
solved to some extent by elaborating the numerical schemes for computation of Pfaf- 
fians [ 12 El in El • The numerically derived^ results for xx {yy, xy, yx) dynamics 
supplemented by the analytical results for zz dynamics [EI permit to work out the 
theory of dynamic properties of spin-| XY chains in the fermionic picture (G.Miiller 
with coworkers) and thus to explain the peculiarities of responses of the spin system 
to small external perturbations [ |Hl E] • 

Let us sketch briefly the linear response theory of the spin chain in the fermionic 
language considering for simplicity the isotropic XY model {J^ = J"^ = J). The zz 
dynamic structure factor for this model which is given by 

/TT 
dmn^^ (1 - n^^_^) 5{uj + A^^ - ^^^-^,) , (2) 
■TT 

1 



A^ = f2 -|- J COS K, 



l + exp(/5A,) 



suggests the following interpretation [El GDI- Consider at first the high-temperature 
limit /3 ^ when — > | and hence ^^^(k,^;) Q becomes independent of VL. Ap- 
plying the infinitesimally small external field (directed along z axis) characterized 
by the wave vector k and frequency uj we observe that the responsive magnetiza- 
tion (directed along z axis) is determined by generation of the two fermions with 
energies A^j and under the restrictions k, = ni — k,2 and u = — A^^ -|- A^j = 

^The numerical approach is not restricted to the uniform chains and can be easily applied to 
the nonuniform chains in which the Hamiltonian parameters vary regularly along the chain with a 
finite period or are random variables with a given probability distribution [|7]. 
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Figure 1. 52z(K,a;) for the isotropic XY chain in a transverse field (J = 1) at 
infinite temperature /3 = (a) and zero temperature /3 — > oo (b: = 0, c: 
n = 0.5, d:n = 0.9). 



2Jsin|sm (ki — l). The "dummy" wave vector ni in Q varies within the region 
—71 ^ ^ TT. As a result, such an experimental probe "measures" a continuum 
of the two-fermion excitations in the k-uj plane. The upper boundary of the two- 
fermion continuum is given hj u = 2| Jsin || and ^'^^(fi;, cu) exhibits divergence along 
this line as it follows from (0). At low temperatures /3 — *■ oo, the Fermi factors in Q 
come into play and Szz{i^,i^) becomes dependent on Q. The additional conditions 
A^i < and A^j > lead to the appearance of the lower boundary (for example, 
uj = |Jsinfi:| if = 0) at which a finite value of S'^^(k,ci;) jumps to zero. The trans- 
verse field (which plays a role of the chemical potential in the fermionic picture) 
effects the lower boundary of the two-fermion continuum in the k-uj plane and the 
redistribution of the values of Szz{i^,i^) in the k-u plane at low temperatures. For 
\Q\ > \J\ Szz{i^,uj) vanishes everywhere in the k-uj plane. The above-said can be 
seen in figure 1 where some typical results illustrating zz dynamics are reported. 

Contrary to the zz dynamics the xx dynamics is more involved: we do not know 
the explicit expression for Sxx{i^,^^) similar to Q. Equation Q arises after compu- 
tation of the average of the product of four Fermi operators 

((l-2c+(t)c„(t)) (l-2c++^c„+^)) 
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that is obviously the two-fermion quantity. The xx dynamic structure factor contains 
the averages hke 

((1 - 2c+(t)ci(t)) ... (1 - 2c+„i(t)c„_i(t)) (c+(t) + c„(t)) 

X (l — 2c^Ci) . . . (l — 2c,^^„_^C„+m-l) (c.^+m + C„+m) ) 

and thus it is a many-fermion quantity. However, the numerical calculations show 
that the two-fermion continuum dominates the low-temperature behaviour of the xx 
dynamic structure factor. Although S^xi''^-,^) is not restricted to the two-fermion 
continuum region in the k-uj plane and has nonzero value above the upper boundary 
of the two-fermion continuum (demonstrating the effects of the many-fermion con- 
tinua [E]) its value outside the two-fermion continuum is rather small. Sxx{f^i^) 
may be described by several washed-out excitation branches following roughly the 
two-fermion continuum boundaries. Further studies are required to clarify why two- 
fermion features rule the many-fermion quantity S^xif^-,^)- For |^| > |<^| the zero- 
temperature XX dynamic structure factor Sxx{.^^^^) shows a single 5-peak along the 
fermion branch uo = A^. In the high-temperature limit the low-temperature struc- 
tures in K-uj plane disappear and Sxx{i^,^) becomes ^-independent in agreement 
with exact calculations for /? = [E]. Alternatively the xx dynamics can be exam- 
ined using a bosonization treatment however, such an analysis is restricted to 
low-energy physics and only a small region in the k-uj plane can be explored by this 
approach. 

The described analysis of the dynamic properties may be extended to the ani- 
sotropic XY interaction (J^' ^ J^) [[S] and the dimerized isotropic XY interaction 
(J^' = — *• J(l — (— ^ 5 ^ 1 is the dimerization parameter) [IHIE]- Appar- 
ently, the exhaustive study of the zz dynamics for the former case [|HI should be still 
supplemented by the corresponding analysis of the xx dynamics whereas the case 
when both anisotropy and dimerization are present requires a separate study. Re- 
cently, the effects of periodic inhomogeneity on the dynamic susceptibility Xzz{,i^,^) 
(but not Xxxif^i^)) have been reported in the paper on dynamics of isotropic XY 
model on one-dimensional superlattices [ E] • The xx dynamic quantities for the 
case of extremely anisotropic XY interaction (J^ = 0), i.e., for the spin-i transverse 
Ising chain, are of interest for interpreting the experimental data on the dynamic 
dielectric permittivity of the quasi-one-dimensional hydrogen-bonded ferroelectric 
compound CSH2PO4 [CnillE]. 

To end up this section, let us note that the dynamic properties of two-dimensional 
quantum spin models can be also explained in terms of the two-fermion continuum, 
however, such a picture may have only an approximate meaning [ El • Another 
interesting question is to contrast the results for dynamic structure factors of spin-| 
and spin-1 chains [flHj. 

3. Dzyaloshinskii-Moriya interaction 

The Dzyaloshinskii-Moriya interaction is often present in the low-dimensional 
quantum magnets (see, for example, a recent paper [EI). It is generally known 
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that the Dzyaloshinskii-Moriya interaction Yin ^ ' i^n ^ ^n+i) being added to the 
Hamihonian (P) does not destroy the rigorous treatment if D = (0, 0, D) [ 20J^. 

The main effect of the Dzyaloshinskii-Moriya interaction is the loss of the sym- 
metry of elementary excitation energies with respect to the change n —>■ —n: 



A«; = D sin K + y \^-\ cos nj + J sin k ^ 

[123 1221 • presence of the Dzyaloshinskii-Moriya interaction some remarkable 

changes in the thermodynamic and dynamic properties of spin chains may occur. 
Consider, for example, the isotropic XY chain. Since the Dzyaloshinskii-Moriya 
interaction for such a chain (even inhomogeneous one) can be eliminated by the 
local rotations in spin space around the z axis 



S"^ cos (pn + Sn Siu 0„ S^, -S^ siu + COS 



n ' 



V9i + . . . + v?n-i, tanv3„ = ^ (3) 



'm 



resulting in a model with isotropic XY interaction \/ Jn + D"^ the zz dynamics 
remains as for a chain without Dzyaloshinskii-Moriya interaction, however, with 
renormalized energy scale J — > V + D"^ ^. In contrast, the xx dynamic quantities 
according to ^ involve 

«(^X+m)j,D = COS (t)n COS (t)n+m{s nit) J ^p^:^^^ 

— COS0„ sin (l)n+m {S nit) Sn+rn)^J^+D^,0 

- sin0„cos0„+^(s;^(tX+^)^AjT^_o 

+ sin0„sin0„+^(s^(t)4+„)^/j2^o (4) 

(for the homogeneous chain 0„ = (n — tany9 = D/ J). In view of ^ the relation 
between the zz and xx dynamics discussed in the previous section (S'^^(fi;,co') at low 
temperatures exhibits the washed-out excitation branches which roughly follow the 
boundaries of the two-fermion continuum which determines 6*22 (/€,u;)) may appear 
to be more intricate [^^. In the next sections 4 and 5 we give further examples of 
how the Dzyaloshinskii-Moriya interaction manifests itself in the properties of spin-| 
XY chains. 



4. Correlated ofr-diagonal and diagonal disorder 

The Jordan- Wigner transformation maps the spin-i isotropic XY chain in a 
transverse field onto the chain of tight-binding spinless fermions with the on-site 

^Let us note that in some cases the Dzyaloshinskii-Moriya interaction can be ehminated by the 
corresponding spin-coordinate transformation [E] (see, for example, equation ^). 

^This is not the case if XY interaction is anisotropic; the zz dynamics of the model with 
extremely anisotropic XY interaction, i.e., the transverse Ising chain with Dzyaloshinskii-Moriya 
interaction, was considered in [\22\. 
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energy Q and hopping / = J/2. If the transverse fields are independent random 
variables (diagonal disorder) each with the Lorentzian probability distribution 



1 



p{VLn) = -— —2 



the resulting fermionic model is the one- dimensional version of the Lloyd model. 
The density of states, 



1 ^ r 

P^^) = ]7 E - ^^)' (•••) = •••/ d^]nP(fin) ■■■(■•■), 

k=l ^ 

for the Lloyd model can be found exactly ["^T. Going far beyond the idea of H.Ni- 
shimori we may consider a spin model with the correlated off-diagonal and diagonal 
Lorentzian disorder which after fermionization reduces to the one-dimensional ver- 
sion of the extended Lloyd model introduced by W.John and J.Schreiber. Namely, 
we consider the isotropic XY model with independent random exchange interactions 
(off-diagonal disorder) given by the Lorentzian probability distribution 

Moreover, we consider the correlated off-diagonal and diagonal disorder assuming 
that the on-site transverse fields in the chain are determined by the surrounding 
exchange interactions according to the relation 

fi„ - ilo = ^ {Jn~\ + Jn- 2Jo) , a is real, \a\ ^ 1. (6) 



Then the density of states p{E) yielding the thermodynamic quantities of the intro- 
duced random spin chain can be calculated exactly [ESlESj- To get p{E) we must 
calculate the diagonal Green functions G^^{E), since 



The set of equations of motion for G^^{E ± ie), e — > +0 can be averaged using 
contour integration in complex planes of random (Lorentzian) variables J„. Using 
the Gershgorin criterion we find the set of equations for the averaged Green func- 
tions which has the same structure as before averaging but possesses translational 
symmetry. As a result we obtain the desired quantities G^^{E) and hence all ther- 
modynamic quantities. 

In figures 2a and 2b we present p{E) in the most interesting region \a\ 1, when 



p{E) becomes not symmetric with respect to the change E — ^Iq — > — {E — ^Iq). Such 
asymmetry immediately yields a nonzero (average) transverse magnetization 



I dEp(E) tanh ^ ^ 
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Figure 2. The density of states (a, b) and the transverse magnetization (c, d) 
of the isotropic XY chain with correlated Lorentzian disorder (Jq = 1, T = 1, 
a = —1 (a, c), a = 1 (b, d)). The dotted curves correspond to the nonrandom 
case (r = 0). 



at low temperatures /3 — oo at zero (average) transverse field f2o = (figures 2c 
and 2d). Let us consider more closely this somewhat unexpected magnetic property 
of the introduced random spin chain. For a certain random realization of the chain 
defined by (0), (jHl) one may expect the same numbers of sites surrounded by stronger 
than Jo exchange interactions as the sites surrounded by weaker than Jq exchange 
interactions. Because of (jHl) for f2o = the transverse fields at the former and at the 
latter sites have the same value but the opposite signs giving as a result ^„ fi„ = 0. 
On the other hand, one may expect that the sites surrounded by strong isotropic XY 
exchange interactions exhibit small z magnetization whereas the sites surrounded 
by weak isotropic XY exchange interactions exhibit large z magnetization (in the 
opposite direction). As a result, the average transverse magnetization has a nonzero 
value. As \a\ increases, a difference in the oppositely directed z magnetizations be- 
comes smaller. Thus, a nonzero at fio = appears owing to the imposed relation 
(0) which expresses the condition of correlated disorder. 
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Some further insight into the origin of the asymmetry of p{E) can be obtained 
after examining the moments of the density of states 



(here H is the Hamiltonian of fermions which represent the spin chain). It is just 
the correlated disorder that yields a nonzero third moment M^^) 7^ 0"^ at fio = 0. For 
not correlated off-diagonal and diagonal disorders one gets M(^) = at fio = 0. The 
moments of the density of states can be calculated for any probability distribution of 
random variables p{Jn) (not necessarily for the Lorentzian probability distribution), 
for example, for the rectangle probability distribution. These calculations explicitly 
demonstrate the cause of the asymmetry appearance in p{E) [EZj. Some other results 
on the effects of correlated disorder can be found in [ HUl I2HI 12^] • 

Finally, let us remark that the considered spin model with correlated Lorentzian 
disorder may be extended by introducing the nonrandom Dzyaloshinskii-Moriya in- 
teraction D. Another extension is to assume the exchange interaction to be nonran- 
dom Jn = J whereas the Dzyaloshinskii-Moriya interactions Dn to be independent 
random Lorentzian variables determining the transverse fields according to (jB)) . Both 
models are related to each other through a certain sequence of rotations of spin axes 
around the z axis (being applied to the Hamiltonian with the exchange interactions 
Dn and the Dzyaloshinskii-Moriya interactions — J„ it gives the Hamiltonian with the 
exchange interactions J„ and the Dzyaloshinskii-Moriya interactions Dn) and hence 
it is sufficient to consider only one of them. Considering, for example, the former 
model one finds that the nonrandom Dzyaloshinskii-Moriya interaction may lead to 
the recovery of the symmetry with respect to the change E — Qq ~ {E ~ Qq) and 
hence to a decrease of the nonzero value of at Qq = 0. Such an effect becomes 
also apparent after calculating the moments of the density of states M^^) and M(^) 
((Zj) for this model. 

To end up this section, let us note that a different random spin-| XY chains 
were rigorously analytically examined by Th.M.Nieuwenhuizen and coauthors [ l3Uj . 

5. Effects of regularly alternating bonds and fields 

5.1. Continued fractions 

The quantum spin chains with regularly alternating Hamiltonian parameters can 
model dimerized (/-merized) chains, one-dimensional superlattices, one-dimensional 

"^Here (...) denotes the random averaging either for the correlated ofF-diagonal and diagonal 
disorder or for the independent off-diagonal (with the probability distribution p(J„)) and diagonal 
(with the probability distribution p(rin) — J dJ„_i J dJnP{Jn-i)p{Jn)5 {p^n — ^^o — f [Jn-i + Jn 
— 2Jo))) disorders. 





n=l 
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decorated chains etc. The case of spin-| XY chains supplemented by the continued 
fraction approach (for other approaches see remains amenable for rigorous 

analysis of the thermodynamic properties if the exchange interaction is isotropic 
(J^ = jy) extremely anisotropic (J^ = 0) [ 33j or if fi„ = [El- The thermo- 
dynamic quantities of the regularly alternating isotropic XY chain in a transverse 
field can be obtained through the density of states p{E) = J2k=i ^{E — A^). The 
thermodynamic quantities of the regularly alternating transverse Ising chain and the 
regularly alternating anisotropic XY chain without transverse field can be obtained 
through the density of states R{E'^) = ^ ^f^^ 6{E^ - Al). 

Let us recall that after fermionization of the isotropic XY chain in a transverse 
field {Jn — Jn ~ '^■^n) oue faces the Hamiltonian which is a bilinear fermion form. 
While making it diagonal one arrives at the set of equations 

In-igk,n-l + (^n — Afc) Qkn + Ingk,n+l = 0, (8) 

here gkn are the coefficients of the linear transformation which make the initial bilin- 
ear fermion form diagonal and are the resulting elementary excitation energies. 
Therefore, introducing the Green functions Qnm{E) according to 

— In-lQn-l,miE) + {E — Qnm{.E) — In'3n+l,m{.E) = 5nm (9) 

one gets the density of states 

1 ^ 

p{E) =T^y2^Gnn{E±ie), 

TCN ^ 

n=l 

For the transverse Ising chain (J^ = 2/„, = 0) instead of (jH)) and Q we have 

nn-lln-l^k,n-l + {^l + " A^) + UJ^^k^n+l = (10) 

and 

— ^ln-lIn-lQn~l,miE'^) + (E'^ — — Qnm{E'^) — VLnInGn+l,m{E'^) = Snm 

respectively, whereas for the anisotropic XY chain without field {J^ = 21^, = 21^, 
fin = 0) instead of (jH)) and © we have 

^24^1'^ ;^,n-2 + (iLl' + If - Al) <^kn + 4^/:+l$fc,n+2 = (12) 

and 

- In-2ln-lSn^2,m{E'^) + [E"^ - In-l " if) Gnm{E^) - I^In+iGn+2,m{E'^) = Snm , 

(13) 

respectively. For the last two models the Green functions Qnm{E'^) yield the density 
of states 

1 ^ 

R{E') = ^ '^GnniE' ± ie), e ^ +0. 

n=l 
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For a general case of the anisotropic XY chain in a transverse field a set of equations 
like (fTUj) or (fT^ is five diagonal banded (but not three diagonal banded as (jH)), 
(fTUI) or (22)) and the next step, i.e., the continued fraction representation for the 
diagonal Green functions is less evident. According to equation Q, (fTT|) or (fT^ 
the diagonal Green functions for all these models can be represented in terms of 
continued fractions. For example, from (jl3p it immediately follows that 



A 



n— 1 " n r 

TV "ijx 2 
-'n-2 -'ra-1 

n ~ TV 2tx 2 

J^2 _ JX 2 _ jy 2 _ n-4, 'n-3 



-3 n—2 jp2 Tx 2 jy 2 

-'n-5 -'n-4 



jy'ijx 2 

TP2 _ Tx 2 _ T-J/ 2 ^„+2 -'n+3 

-'^ -"n+l -'n+2 7^2 2 fy 2 

-'^ ~-'n+3 ~-'n+4 ~ 

If now the Hamiltonian parameters are periodic with any finite period p the contin- 
ued fractions AJ in (fT^ become periodic and can be evaluated by solving square 
equations. As a result, we get the exact expressions for (or Qnn{E) for the 

isotropic XY chain) and hence for all thermodynamic quantities of the regularly 
alternating spin chains in question. 

Let us remark, that the thermodynamic quantities of the anisotropic XY chain 
in a transverse field of period 2 were obtained in [ in a different manner on 
a quite general background. This scheme, however, becomes cumbersome if the 
period of nonuniformity increases. Further development of a general approach to 
the study of thermodynamic quantities and spin correlation functions has been re- 
ported by L.L.Gongalves with coworkers in connection with spin-i XY models on 
one-dimensional superlattices. 

5.2. Magnetization processes 

The most spectacular manifestation of the effects of regular alternation can 
be seen in the magnetization processes at low temperatures. Consider at first the 
isotropic XY chain. Owing to regularly alternating parameters of the Hamiltonian 
the fermion band splits into several subbands the number of which does not exceed 
the period of regular inhomogeneity. This circumstance immediately suggests that 
the zero-temperature dependence 

1 f j3E 

rriz = — / dEp{E) tanh 

2, J 2 

on Q (we have assumed that f2„ = f2 + AQn) consists of horizontal parts (plateaus) 
which appear when E = remains in between subbands (given by p{E)) with 
changing of Q separated by the parts of varying which appear when E = 
remains within subbands with Q changing. Clearly, the number of plateaus does 
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not exceed the period of nonuniformity. Moreover, their heights —s ^ ^ s 
are in agreement with the condition: p^s — niz) = integer with s = ^ [|36j. This 
conjecture was suggested on the model-independent background for a general spin- 
s chain with axial symmetry in a uniform magnetic field using the Lieb-Schultz- 
Mattis theorem and the bosonization arguments. Obviously, spin-^ XY chain is 
an easy case (represented by noninteracting fermions) contrary to the chains with 
Heisenberg exchange interaction or higher spins s > |. On the other hand, the 
magnetization curve for spin-i XY chain can be obtained explicitly. Moreover, the 
elaborated scheme also permits to obtain the local (on-site) magnetizations. The on- 
site magnetizations exhibit plateaus which begin and end up at the same values of Q 
as for the total magnetization, however, the plateaus heights (i.e., the values of the 
on-site magnetizations) are not universal quantities no longer obeying the above- 
mentioned condition and strongly depend on a concrete set of the Hamiltonian 
parameters. Moreover, a sequence of sites rii, 77-2, rip satisfying the inequalities 
nT'zni ^ '^2n2 ^ • • • ^ ^zn dcpcuds ou the value of the applied field fl. The step-like 
behaviour of the zero-temperature dependence (^zn) vs. Q is accompanied by 
the corresponding zero-temperature behaviour Xz = drriz/dQ {Xzn = dmzn/dfl) vs. 
fl. The zero-temperature dependence Xz on Q behaves like —p{E = 0) vs. Q thus 
reproducing the density of states p{E). 

In contrast to the isotropic XY chain, a regular alternation of the parameters 
of the transverse Ising chain Hamiltonian does not lead to plateaus in the zero- 
temperature dependence on Q. The difference arises owing to the anisotropy 
of XY interaction: if 7^ one finds that [J2n ^] 0- Thus, even if the 
spin system remains in the same ground state |0) with varying Q the magnetization 
(0|^^s^|0), nevertheless, varies with varying fl. As a result the dependence 
vs. Q does not exhibit horizontal parts. On the other hand, in the vicinity of some 
fields Q* (critical fields) the magnetization may behave as ~ (^7 — Q*) ln\Q — Q*\ 
and hence the susceptibility as Xz ~ ln|i7 — For the uniform chain there are 
two critical fields Q* = ±|/|. If the bonds become regularly alternating only the 
values of f2* change. For example, f2* = ±^/\LlI^ for p = 2, Q* = ±-^1/1/2/3! for 
p = 3 etc. If the fields become regularly alternating not only the values of the critical 
fields vary but the number of the critical fields may change. For example, assuming 
^1,2 = ± Afl, AQ ^ one has either two critical fields = ±a/ AQ^ + Ihhl 
if An < ^/\hh\ or four critical fields fi^ = ±^/A^F±~\hh\ if Afi > v^T^^. 
The transverse field is a controlling parameter of the second-order quantum phase 
transition in the spin-| Ising chain in a transverse field [EZl- Thus, a number of the 
quantum phase transition points governed by fl in this model may increase due to a 
regular alternation of the transverse fields presuming the strength of inhomogeneity 
is sufficiently strong. For example, for p = 2 there may be either two or four quantum 
phase transition points, for p = 3 there may be either two, or four, or six quantum 
phase transition points etc. The critical behaviour remains unchanged and is just 
like for the temperature driven phase transition in the square-lattice Ising model. 
For example, the zero-temperature dependence Xz vs. Q always exhibits logarithmic 
singularities the number of which depends on a concrete set of the Hamiltonian 
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parameters. It seems interesting to trace the changes in magnetization processes as 
anisotropy of the exchange interaction varies between the isotropic and extremely 
anisotropic hmits. 

It is instructive to compare the zero- 
temperature magnetization processes in the 
quantum and classical XY chains^. A classical 
chain consists of arrows (vectors) s = {s,6,(j)), 
s = l, 0^6*^7?, O<0<27r which in- 
teract with each other and an external (trans- 
verse) field. The classical isotropic XY chain is 
described by the Hamiltonian 



H = ris cos On 

n 

+ Js^ sin 9n sin 9n+i cos i 



(15) 



(compare with (^). The Hamiltonian of 
the classical transverse Ising chain contains 
cos cos instead of cos (0„ — 0n+i) in the 
interspin interaction terms in (jl5p . Consider fur- 
ther a chain of period 2 with Qi^2 = ± Afi, 
AQ ^ 0, Ji,2 = J- The ground-state energy 
ansatz for both the isotropic XY and transverse 
Ising chains reads 

N 

^0(^1,^2) = -s{{n + An) cosOi 

+ {n - An) 008 02) - N\J\s^ sin 01 sin 02 (16) 

and 01, 02 are determined from the equations 

(fi + Afi) sin6'i + 2| J|scos6'i sin6'2 = 0, 
(fi- Afi)sin^2 + 2|J|ssin0icos02 = (17) 

to provide a minimum of £'0(^1,^2) (|16|). The 
on-site magnetizations have the component ar- 
bitrary directed in the xy plane with = 
ssin^n (n = 1,2) for the isotropic XY case or 
the component directed along the x axis with 
I'fTT'xnl = ssin0n for the transverse Ising case 
and the z component given by = scos6'„ in 
both cases. The components of the total mag- 
netization per site are then as follows = 




Figure 3. The zero-temperature 
magnetization curves for the 
isotropic XY (a), transverse 
Ising (b) and classical (c) chains 
of period 2 (J„ = 1, r2i^2 = 
± An, An = (sohd curves), 
AO, = 0.5 (long-dashed curves), 
AO, = 1 (short-dashed curves), 
AQ = 1.5 (dotted curves)). 



^Rigorous calculations of the initial static susceptibility tensor for another nonuniform classical 
chain, i.e., the spin-i Ising chain, can be found in [I38|. 
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I (m_Li + m_L2) or m^. = i {rrixi + 171^2) (if ^ < 0) and = | (m^i + m^g). Equa- 
tions ()17p may be solved analytically yielding the ground-state energy and the on- 
site magnetizations m±i (mj^i), m±2 ("^0:2)) nT'zi, ^z2- In figure 3 we contrast the 
zero-temperature magnetization curves for the isotropic XY, transverse Ising and 
classical chains of period 2. Note, that the classical isotropic XY chain similarly 
to its quantum counterpart may exhibit plateaus in the dependence vs. Q. Ev- 
idently, it would be interesting to examine a quantum-to-classical crossover in the 
magnetization curves. Probably a real-space renormalization-group method aimed 
on the study of quantum fields on a lattice and applied in [EH] to the spin-i trans- 
verse Ising chain may be used to analyse the spin-1 chain, that is, an intermediate 
case between the quantum (s = |) and the classical (s — > 00) cases. 

Exact results for the magnetization curves of spin-i XY chains may be employed 
to test the strong-coupling approach that is a widely used approximate approach in 
the theory of low-dimensional spin systems consisting of periodically repeating units 
(spin chains with a periodic modulation of the intersite interactions or spin ladders) 

5.3. Spin-Peierls instability 

Since the paper of P.Pincus [E] we know that the spin-i isotropic XY chain is 
unstable with respect to a lattice dimerization due to the spin-Peierls mechanism^. 
Really, the ground-state energy per site of the dimerized spin chain (i.e., with J„ = 
J (1 — (— is given by 

eo(5) = d^^l-(l-52)sin2^ _ \n\ - ^) , (18) 

r 0, if \j\ ^ 

^ = I arcsin^i^p^, if 5\J\ ^ \n\ < \ J\, 
[ f, if 1^1 < 5\J\. 

At n = the ground-state energy per site (HH)) (eo(5) ^ eo(0) + ^^^In^, 5 < 1) 
decreases rapidly enough in comparison with the increase of the elastic energy per 
site aS"^ as 6 increases to provide a minimum of the total energy per site eo(5) + 
at a nonzero value of the dimerization parameter 6*. The external field may destroy 
dimerization: if Q exceeds a certain value the dimerized phase cannot survive and the 
uniform phase becomes favourable. The phase diagram shown in figure 4a specifies 
the region of stability/metastability of the dimerized phase at zero temperature. 

It is generally known (see, for example, a review on CuGeOs [ESI) that the in- 
crease of external field leads to a transition from the dimerized phase to the incom- 
mensurate phase rather than to the uniform phase that contradicts to what is seen 
in figure 4. Obviously, the incommensurate phase cannot appear within the frames 
of the adopted ansatz for lattice distortion 6162S162 . . . , 61 + 62 = 0. Therefore, we 

^This conclusion was obtained within the adiabatic approximation; the corresponding treatment 
within the nonadiabatic approximation is more sophisticated [ I42| . 



741 



O.Derzhko 



should assume J„ = J (1 + 5„) and exam- 
ine the total energy Eq ({5n}) + ^„ ctS^ for 
different lattice distortions We intro- 

duce a trial distortion of the form 



5n = —S COS ( — n 

V P 



(19) 



where p is the period of modulation and 
analyse numerically (A^ = 1000) the total 
energy to reveal which spin-Peierls phase is 
realized in the presence of external field (for 
details see [EI). Comparing the behaviour 
of the total energy for p = 2 and p = 1.9, 
p = 2.1 as Q increases we conclude that 
a long-period structure does arise if Q ex- 
ceeds a certain value. The dimerized phase 
transforms into a long-period phase rather 
than into the uniform phase while the field 
increases. The important conclusion of fur- 
ther computations using (fT^ is that the 
dimerized phase persists up to a certain 
characteristic field. Further, in moderate 
fields the lattice parameterized by (fT^ may 
exhibit short-period phases, for example, 
the trimerized phase with p = 3 [ ES] (an- 
other possible lattice distortion which pre- 
serves the chain length is 6i = —62)- How- 
ever, a behaviour of the trimerized phase is 
essentially different in comparison with that 
of the dimerized phase: for any small devia- 
tion of the field from the value at which the 
trimerized phase occurs there exists such a 
long-periodic structure for which the lat- 
tice distortion (fT^ gives smaller total en- 
ergy than for p = 3. Thus, contrary to 
the dimerized phase, the trimerized phase 
does not persist with the field varying and 
it continuously transforms into a certain 
long-period phase with the field varying. In 
strong fields the uniform lattice can be ex- 
pected. Clearly, such a study is restricted 
to the adopted ansatz (fT^ (although any 
other (but not all) distortion pattern may 
be assumed) and therefore we can say for 
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Figure 4. The phase diagram of the 
isotropic XY chain (J = 1) in the 
plane lattice stiffness a - transverse 
field J7 which indicates the regions of 
stability/metastability of the dimer- 
ized and uniform phases at zero tem- 
perature. A (C) - only dimerized 
(uniform) phase occurs; Bi (B2) both 
phases are possible but the dimerized 
(uniform) is the favourable one. The 
phase diagram also illustrates the 
possible effects of the Dzyaloshinskii- 
Moriya interaction (a: D = 0, h: D = 
0.5, k = l,c:D = 0.5, k = 0). 
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sure which lattice distortion is not reahzed rather than to point out which lattice 
distortion should occur. 

It is interesting to discuss the effect of the Dzyaloshinskii-Moriya interaction 
on the spin-Peierls dimerization. For this purpose we should find the ground-state 
energy of the chain of period 2 with a sequence of parameters 



Moreover, Ji^2 = ^(1 ± ^) and Di^2 = -0(1 ± k6). Putting k = one has a chain in 
which D does not depend on the lattice distortion, whereas for k = 1 the dependence 
of D on the lattice distortion is the same as that for the isotropic XY exchange 
interaction. In the limit of small 5 -C 1 (valid for lattices with large values of a 
and corresponding to the experimental situation) the analysis becomes extremely 
simple. It is convenient to introduce the parameter 



Consider, for example, the case f2 = 0. After a simple rescaling of variables one 
arrives at the equations considered in [HTj and as a result 



dimerization parameter 6*, whereas for k = nonzero D leads to the decrease 
of 6*. Thus, the Dzyaloshinskii-Moriya interaction may act either in favour of the 
dimerization or against it. The actual result of its effect depends on the dependence 
of the Dzyaloshinskii-Moriya interaction on the amplitude of the lattice distortion 
in comparison with the corresponding dependence of the isotropic XY exchange 
interaction. Comparing figure 4a and figures 4b, 4c one can see the described effects 
of the Dzyaloshinskii-Moriya interaction on the spin-Peierls dimerization. 

Finally, let us make some remarks on the exchange interaction anisotropy effects 
on the spin-Peierls dimerization. Consider the anisotropic XY chain (without field). 
In the isotropic limit it exhibits the spin-Peierls dimerized phase, whereas in the 
extremely anisotropic limit (Ising chain) the ground-state energy does not depend 
on 6 and the spin-Peierls dimerization cannot occur. Therefore, it is attracting to 
follow how anisotropy of the exchange interaction destroys the dimerized phase 
and to find, in particular, a critical value of anisotropy above which the dimerized 
phase cannot appear. The consideration becomes even more complicated in the 
presence of an external (transverse) field. Note, however, that such a study could be 
performed many years ago on the basis of the results for the anisotropic XY chain 
in a transverse field of period 2 reported in [ I35j . 

Another issue which deserves to be studied in detail is related to a spin-Peierls 
instability of the spin-i transverse Ising chain [ISj- By direct inspection one can 
make sure that although the ground-state energy of the dimerized Ising chain in a 



njiDiQJ2D2^JlDiQJ2D2 ... . 



J2 + ■ 




We immediately conclude that for k = 1 nonzero D leads to the increase of the 
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constant transverse field Q decreases as S increases, however, not sufficiently rapidly 
(for any Q) to get a gain in the total energy. As a result no spin-Peierls dimerization 
should be anticipated. On the other hand, the transverse Ising chain is unitary 
equivalent to the anisotropic XY chain without field. By comparison of the diagonal 
Green functions (see equation (fT^ in the present paper and equation (8) in [ EHl 
J. Phys. A]) it can be found that the Helmholtz free energy of the anisotropic XY 
chain without field of N sites J^J^J^_,_^J^_,_]^ ... is the sum of the Helmholtz free 
energies of two transverse Ising chains flnln^n+iln+i • • • of N/2 sites, namely, 

TX jy TX jy 

and 

jy TX TV TX 

■ ■ ■ -'n-l-'n-'n+l-'n+2 • • • • 

Consider at ffist the following isotropic XY chain = = / (I — (— which 
exhibits the spin-Peierls dimerization. In view of the mentioned correspondence the 
introduced model is thermodynamically equivalent to two uniform transverse Ising 
chains (both with the transverse field /(I + S) (or /(I — 6)) and the Ising exchange 
interaction — S) (or + S))) and thus we conclude that the critical Ising chain 
[Q = I) is unstable with respect to a uniform extending (or shortening) accompa- 
nied by the corresponding increase (decrease) of the transverse field. Further, it is 
known that the quadrimerized isotropic XY chain without field may be energetically 
favourable in comparison with the uniform chain (although yielding a smaller gain 
in the total energy than the dimerized chain). Consider, therefore, such a chain of 
period 4 with J„ = /'(I + 6), h+i = /"(I + 5), J„+2 = /'(I - 5), 4+3 = /"(I - ^) 
which is unstable with respect to a lattice distortion characterizing by nonzero value 
of 5 for certain values of J', J". This model exhibits the same thermodynamics as 
two identical transverse Ising chains of period 2 f2„/„f2„+i/„_|_ifi„/„i7„_|_iJ„+i . . . with 
fi„ = J"(l-5), Vtn+i = /"(1+5) and J„ = /'(1+5), In+i = These arguments 

indicate a possibility of the spin-Peierls bond dimerization of the spin-| transverse 
Ising chain accompanied by the coherent modulation of the on-site transverse fields. 



6. Skipped items and summary 

Needless to say that many important contributions have appeared out of the 
scope of this brief survey. Let us discuss in a telegraph-style manner some of them. 
Numerous works were devoted to the analysis of the properties of spin-| XY chains 
with aperiodic or random Hamiltonian parameters [,46,; (as, for example, an extensive 
real-space renormalization-group treatment of the random transverse Ising chain by 
D.S.Fisher). 

One of the interesting subjects in the theory of magnetic materials is the magnetic 
relaxation in spin systems and, in particular, the impurity spin relaxation in spin 
systems. Spin-i XY chains provide a possibility to rigorously study the relaxation 
phenomena in one- dimensional spin models containing impurities [ HZl Ej • A single 
impurity may be assumed in the sense that the interaction with its neighbours is 
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different in strength [J' ^ J). The impurity spin may be located either at the 
boundary of the system (Ji = J', J2 = J3 = . . . = J) or in the bulk (e.g., Jn_^ = 
Jn_ = J' and Jn = J for all other n). After the Jordan- Wigner transformation one 
faces a chain of tight-binding fermions with the impurity site (in the sense that the 
hopping amplitude (s) surrounding this site is (are) different in value) the energy 
spectrum of which is well-known (if J' < the elementary excitation energies form 
a band whereas if J' > two states emerge from the band; moreover, 
[Jc = J) for the boundary (bulk) impurity). Nevertheless, the time dependence of 
the equilibrium autocorrelation functions (s"(t)s") is not obvious since they are 
two-fermion (zz) or, generally speaking, many-fermion (xx) quantities. The time- 
dependent autocorrelation functions (s"(t)s") for the impurity spin or for the spins 
in its vicinity may exhibit new types of asymptotic behaviour depending on the 
relation between J and J' and temperature. The impurity relaxation may become 
even more complicated in the presence of an external (transverse) field. 

Spin- 1 XY chains have been used to study the nonequilibrium properties of 
quantum systems [IIHI- Z.Racz with coworkers suggested to consider a nonequilib- 
rium system imposing a current on a system and investigating the steady states. 
Alternatively, we can gain an understanding of the nonequlibrium properties exam- 
ining the dynamics of an initial state (for example, a kink or droplet configuration 
of z on-site magnetizations). 

Recently a study of the transport properties (for example, of the thermal conduc- 
tivity) of low-dimensional spin systems which are significantly determined by mag- 
netic excitations have attracted much interest [llHl- A simple case of the isotropic 
XY chain emerges in such studies as a milestone providing reference results for more 
sophisticated models. 

Finally, some exotic applications of spin-| XY chains have appeared recently. For 
example, in connection with quantum information processing the numerical/analytical 
computations of entanglement in spin-| XY models on one-dimensional lattices of 
small/infinite number of sites were carried out [ED]- A study of the correlation 
function which is called the emptiness formation probability (i.e., the probability 
of the formation of a ferromagnetic string in the antiferromagnetic ground state) 
yields interesting links between statistical mechanics and number theory. The case 
of the isotropic XY chain provides a valuable background for calculation of such 
correlation functions [i51j. 

We hope that the present review shows that spin-i XY chains still contain quite 
unexplored properties which deserve to be discussed. A serious advantage of this 
type of models is a possibility to perform statistical mechanical calculations either 
rigorously analytically or exactly numerically (considering as long chains as required 
to obtain results which pertain to the thermodynamic limit). On the other hand, it 
is always desirable to clarify afterwards the relation of the obtained results to more 
realistic models (with the Heisenberg interaction, interchain interaction, s > |). 
Many interesting questions in the theory of these quantum spin chains still remain 
open and call for new efforts. 
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